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ABSTRACT 



We model the time- variable absorption of Fe II, Fe III, Si II, C n and Cr n detected in UVES spectra of GRB 080310, with 
the afterglow radiation exciting and ionizing the interstellar medium in the host galaxy at a redshift of z = 2.42743. To 
estimate the rest-frame afterglow brightness as a function of time, we use a combination of the optical VRI photometry 
obtained by the RAPTOR-T telescope array - which are presented in this paper - and Swift's X-Ray Telescope 
observations. Excitation alone, which has been successfully applied for a handful of other GRBs, fails to describe the 
observed column-density evolution in the case of GRB 080310. Inclusion of ionization is required to explain the column- 
density decrease of all observed Fell levels (including the ground state 6 Dg/ 2 ) and increase of the Fein 7 S 3 level. The 
large population of ions in this latter level (up to 10% of all Fein) can only be explained through ionization of Fell, 
whereby a large fraction of the ionized Fell ions - we calculate 31% using the Flexible Atomic (FAC) and Cowan codes 

- initially populate the 7 S3 level of Fein rather than the ground state. This channel for producing a significant Fein 
7 S3 level population may be relevant for other objects in which absorption lines from this level - the UV34 triplet 

- are observed, such as BAL quasars and r\ Carinae. This provides conclusive evidence for time-variable ionization 
in the circumburst medium, which to date has not been convincingly detected. However, the best-fit distance of the 
neutral absorbing cloud to the GRB is 200-400 pc, i.e. similar to GRB-absorber distance estimates for GRBs without 
any evidence for ionization. We find that the presence of time- varying ionization in GRB 080310 is likely due to a 
combination of the super-solar iron abundance ([Fe/H]=+0.2) and the low Hi column density (log iV(Hl) = 18.7) in 
the host of GRB 080310. Finally, the modelling provides indications for the presence of an additional cloud at 10-50 pc 
from the GRB with log JV(Hi) ~ 19-20 before the burst and which became fully ionized by the radiation released 
during the first few tens of minutes after the GRB. 
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1. Introduction 



X 

S-H 



predicted, such as th e gradual ionization of H i and Mgii 



Gamma-ray burst (GRB) afterglows can be dete cted at 
nearly any wavelength up to ver y high redshifts ( Tanvir 
|et al. |2009| |Cucchiara et al.||2011|) and are associat ed with 



( Perna & Loeb 1998 Perna & Lazzati 2002[), the excita- 



tEe death s of m assive stars (for a recent review, see Hjorth 



tion and diss ociation of H2 molecules (|Drainc 2000 Draine 



fe Hao 
Draine 



2000 



2002 ) the destruction of d ust grains ( |Waxman fe" 



Fruchter et al.||2001|) and the accompanying 



decrease in extinction and release of metals into the gas 



fe Bloom [201 \\ \ they are therefore considered promising phase (|Perna fe Lazzati 2002 Perna et al.|2003[). Apart from 



probes of star formation at high redshift (e.g. Lamb & 



Reichart 2000 1. However, in order to interpret the wealth of 
information on the interstellar medium (ISM) of GRB host 
galaxies gathered from GRB afterglow sp ectroscopy (e.g 



Prochaska et al. 2007 Fynbo et al. 2009), it is important 



to understand in what way, and up to which distance, a 
GRB explosion is affecting its host. 

Several possible effects due to the brief but extremely 
powerful radiation of a GRB and its afterglow have been 



the detection of excited H 2 molecules ( Shelter et al.||2009 ), 
none of these effects have been convincingly detected. For 
dust destruction, this may be explained by the time scale 
being too short (tens of seconds) for present observations 
to allow a firm detection. 



One effect that was not predicted, but which has 
now been unambiguously observed in several GRBs, is 
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absorption-line variabi lity of fine-structure l inej^l of ions 



such a s Fen and Nili (Vreeswijk et al. 2007 
2009a). This variability has been shown to 



D'Elia et al 
De due to the 

afterglow ultraviolet (UV) photons exciting the neutral ab- 
sorbers in the ISM at distances of a hundred parsec up to 
well over a kiloparsec from the GRB (|Prochaska et al.|2 006 



Vreeswijk et al.|2007l |D'Elia et al.|2009a| ). These distances 



are consistent with lower limit estimates for the neutral 
gas (> 50-100 pc) based on the presenc e of Mgl in the 
afterglow spectra (Prochaska et al. 20061. They are also 



in agreement with hydrodynamic calculations of the size 
of the pre-GKB ionization bubble that is being created by 
the GRB progenitor star and its likely cluster companions 
( |Whalen et al1[2008| . Such a scenario, in which the imme- 
diate environment is already mostly ionized by the time 
that the GRB occurs, can also explain the difference be- 
tween the equivalent hydrogen column density measured 
from the soft X-rays, N(H), and the neutral hydrogen col- 
umn density, iV(Hi), inferred from Lyq absorption in the 
optical/UV spectra ( Wat son et al.||2007| |Campana et al. 
20101 |Schady et al.|[20TT| 



In a companion paper (De Cia et al. 2012 hereafter 
referred to as Paper I), we report in detail on our time- 
resolved high-resolution spectroscopic observations of the 
GRB 080310 afterglow with the Ultraviolet and Visual 
Echelle Spectrograph (UVES), mounted on the Kueyen unit 
of ESO's Very Large Telescope (VLT). This sightline dis- 
plays an unusually low H I column density at the GRB 
redshift (z = 2.42743), log AT(Hi) = 18.7, with an ex- 
treme iron and chromium overabundance: [Fe/H] = +0.2 
and [Cr/H] = +0.7 - these estimates include a correction 
for ionization effects. The values for the carbon, oxygen 
and silicon abundances are instead rather typical for GRB 
sightlincs. Another outstanding feature of the GRB 080310 
UVES spectra reported in Paper I is the unique detection 
of the Fern UV34 triplet at 1895 A, 1914 A, and 1926 A, 
never seen before in a GRB sightline. This, combined with 
the simultaneous decrease of the column density population 
of all levels of Fen, including the ground state, is suggestive 
of on-going ionization at the time the UVES spectra were 
being secured. 

In this follow-up paper we study this hypothesis in de- 
tail, by modelling the column densities of Hi, Fen, Fein, 
Si II, C II and Cr n observed in Paper I as a function of time, 
incorporating - for the first time - both photo-excitation 
and -ionization in a consistent manner. An important input 
parameter for our calculations is the afterglow brightness 
as a function of time, which we estimate by combining the 
observed optical and X-ray fluxes. The latter are derived 
from observations by t he Swift X-Ray Tel escope, which are 
publicly available (see |Evans et al. 2009). For the optical 
we use the clear-filter and VR1 light curves as measured by 
the RAPTOR-T array, which started imaging th e field as 
early as 32 seconds after the GRB trigger time (see Wozniak 
et al. |[20CJ8 ); these data are also presented in this paper. 

This paper is organized as follows. We first describe 
the RAPTOR-T measurements and present the broad-band 
VRI and clear-filter light curves in Sect. [2] In Sect. [3] we 



1 The interaction of the total electron spin and the total elec- 
tron angular momentum causes a fine-structure splitting of the 
atom levels, and the transitions with the lower energy levels cor- 
resp onding to these excite d levels are called fine-structure lines 
(see |Bahcall fc WbTf|l968 1 . 



describe the implementation of the excitation and ioniza- 
tion processes in our modelling code, where the reader is 
referred to Appendix |A"| for the details. We present the re- 
sults of our model fits to the GRB 080310 ionic column 
densities published in Paper I in Sect. [4] These results are 
discussed in Sect. [5] and we briefly summarize our findings 
in Sect. [6] 

2. The RAPTOR-T light curves 

GRB 080310 triggered the B urst Alert Telescope (BA T) 
onboard the Swift satellite (Cummings et al. 2008) at 
08:37:58.65 UT on March 10, 2008. The RAPTORTTtele- 
scope array began observing the BAT position within 10.7 
seconds after receiving the GCN alert, i.e. 32.4 seconds af- 
ter trigger time. RAPTOR-T consists of four co-aligned 
0.4-m telescopes on a single fast-slewing mount and pro- 
vides simultaneous images in four photometric bands (V, 
R, I, and clear). The system - owned and operated by the 
Los Alamos National Laboratory (LANL) - is located at 
the Fenton Hill Observatory at an altitude of 2500 m in the 
Jemez Mountains of northern New Mexico. The RAPTOR- 
T response sequence consists of 9, 20, and 170 exposures 
lasting, correspondingly, 5, 10, and 30 seconds each, and 
separated by 5-second intervals for readout. Approximately 
25% of individual frames were rejected due to intermittent 
glitches in telescope tracking. 

Aperture photometry was per formed on co-added im- 
ages using the Sextractor package ( |Bertin fc Arnouts|l 996) 
with object coordinates fixed at the values measured on the 
reference image where the GRB is detected at a high signal- 
to- noise ratio (S/N) in all four channels. Instrumental light 
curves were then transformed to standard Johnson magni- 
tudes using SPSS photo metry of stars in the vicinity of the 
burst (Cool et al. 20081 and equations of Lupton (2005^] 
The results are listed in Table [T] and plotted in Fig. [T] 



2 http:/ /www. sdss.org/dr5/algorithms/sdssUBVRI 
Transform. html#Lupton2005 
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Fig. 1. Light curves in V (bottom), R, I (top) and the clear 
filter as recorded by the RAPTOR-T telescope. The magni- 
tudes have not been corrected for the Galactic foreground 
extinction. The four epochs at which the UVES spectra 
were taken are indicated with the grey vertical bands. 
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Table 1. Log of RAPTOR- T observations. 
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Notes. (<l) Effective time at measurement midexposurc, start and end, since the Swift BAT trigger on 10 March 2008 at 08:37:58.65 UT. (6) The 
magnitudes have not been corrected for the Galactic foreground extinction. ^ The limiting magnitudes are 3a. 



No optical emission was detected during the first two 
minutes after the burst, down to a limiting magnitude of 
R ~ 18.9 (3<j). The GRB is clearly detected in all sub- 
sequent co-adds starting at 133 seconds after the trigger. 
Following a rapid increase in brightness to a peak value 
at R ~ 16.6 mag, the optical emission fluctuates by a few 
tenths of a magnitude, and after ~ 30 minutes begins a 
slow decline. While the BAT light curve still shows a de- 
tectable gamma-ray e mission between 150 and 320 seconds 
after the trigger (see Littlejohns et al. 2012), the optical 
emission over this time interval appears uncorrelated with 
the gamma rays. The optical light curves from RAPTOR- 
T show no significant color evolution. We determined the 
spectral slope (3 (with F v a v* 3 ) as a function of time, by 
fitting the VRI magnitudes - after correcting the m for the 
Galactic fo reground extinction of Eb-v — 0.045 (Schlegel 
et al. ||1998 1 - at each epoch. The resulting slope values do 



not snow any trend in time, and cluster around the value 
(3 = —1.0, with a standard deviation of 0.4 and an error in 
the mean of 0.07. 

The RAPTOR-T measu rements are gene r ally c onsis 



tent with those reported in Littlejohns et al.| ( 2012[ ). The 
RAPTOR V- and i?-band magnitude limits (3cr), at a mid- 
exposure time of 84 seconds after the burst, correspond to 
Fy < 136 [xJy and Fr < 78 ^Jy, respectively. This F-band 
limit is in agreement with the Swift y-band measurement 
of Fy = 247 ± 140 /iJy at the same epoch, but our i?-band 
limit is well below (almost a fact or of three) the prompt - 
emission fit featured in figure 10 of Littlejohns et al. ( 2012 1 . 



Instead, it is fully consistent with the alternative afterglow 
fit shown in their figure 8. 



3. Modelling the absorption-line variability 

In Paper I, we presented Voigt-profile fits to the aborption 
lines detected in the GRB 080310 spectra, using four dif- 
ferent velocity components. Since these are very close in 
velocity (< 60 km s^ 1 ), it is difficult to ascertain that they 
are indeed correctly separated in the Voigt profile fit, even 
though a strong case can be made that components "b" 
and "c+d" probably are. Moreover, the decomposition is 
not unique, as additional components may be present that 
are hidden in the profile. An added complication is that it 
is unclear which fraction of the H I column density belongs 
to which velocity component, which is important for the 
modelling when ionizing radiation is included. Inspection of 
the Fe n and Fe in column density evolution of the separate 
components indicates a generally similar behaviour, which 
suggests that the components are at a comparable distance. 
Preliminary modelling of the separate velocity components 
"b" and "c+d" indeed results in distances that are the same 
within the error margins. For these reasons, we have focused 
on modelling the total column densities (listed in the last 
column of table 3 in Paper I) rather than those of the in- 
dividual components. 
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3.1. Photo-excitation 

Since the absorption-line variability observed in a handful 
of GRBs can be generally well described by excitation of 
the host-galaxy ISM by afterglow UV photons, we first set 
out to model the GRB 080310 observed column density evo- 
lution as reported in Paper I with photo-excitation alone. 
The excitation fit ting procedure applied here is similar to 
that described in Vreeswijk et al.| ( |2007 ), and which was 
also applied b y Ledoux et al. ( 2009 ) and independently by 



D'Elia et al.| | |2009a|b[ |20*Tof , but with two major differ- 
ences. First, we includ e the correct excitation flux (see the 
erratum published by Vreeswijk et al. 2011), resulting in 
a distance decrease of y/An with respect to the old exci- 
tation calculation. And second, instead of calculating the 
excitation at line centre only, we now effectively integrate 
over the full line profile to obtain the actual flux that is 
entering a particular layer. This second change leads to a 
modest increase in the distance estimate of about 10%. The 
details of the excitation implementation are described in 
Appendix [X] 

For the Fell ion, we use the transition probabilities of 



the 371-level model atom as collected by [Verner (19991, 
including the 63 lower even-parity and 227 higher odd- 
parity levels. Note that the transitions between even-parity 
levels are forbidden and have low transition probabilities, 
while the transitions between even- and odd-parity levels 
are electric dipole, i.e. allowed transitions. The observed 
resonance and fine-structure lines observed in the spectra 
of GRB 080310 and other GRBs correspond to this latter 
group. This 371-level model Fe II atom is sup plemented with 
trans itions taken from the Kurucz data base ( Kurucz & Bell 
|l995[ pj For Fein, we apply two different model atoms. One 
is calculated by Raassen & Uylinga^J from which we include 
59 even-parity, and 214 odd-parity levels. The alternative 
model, which we find to provide a slightly better fit, com- 
bines the A-values for the forbidden transitions between 



the lowest 34 levels calculated by Bautista et al. (20101 



with the allowed transitions from Deb & Hibbert ( 2009p . 



In Paper I we also present measurements and upper lim 
its of the excited-level column densities of Sin and Cn, 
which we include in our fit as well. The t ransition prob- 
abilities of both these ions are taken from Morton ( 2003 ) 
if present therein, and otherwise from the NIST Atomic 
Spectra Databasaj (for the Si n and C II transition probabil- 
ities and their refere nces, see Kelleher fe Podobedova| 2008 ; 



Wiese & Fuhr 
have also calcu 



2007}. We note that |Bautista et aL[ ( 2009| 
ated the A-values of several Si n transitions 



we include the ground level, 2 P°/ 2 
fine-structure level 2 P^ 2 



we find that using those values instead of the NIST ones 
leads to a very similar amount of Si n excitation. For Si n 

and its corresponding 

and 19 higher even-parity levels. 
For C II, we include the lowest two levels ( 2 P°/ 2 an d 2 P^/ 2 )^ 
and 28 higher levels. For both these ions, the energy level 
of the third odd-parity level is larger than the lower even- 
parity levels, so the amount of atoms populating these other 
Sill and Cn odd-parity levels is expected to be negligible. 
Finally, for Cm, we include the lowest 74 even-parity lev- 
els and nearly 400 odd-parity lev els, adopting t he A-values 
of the forbidden transitions from Quinet| ( 1997[ ), and those 



3 http:/ /kurucz. harvard.edu/atoms. html 

4 http:/ /www. science. uva.nl/research/atom/levels/orth/iron 



of the allowed transitions from Kurucz & Bell (1995). For 
all ions, we have made sure that the oscillator strengths 
(or equivalently, the transition probabilities) of the rele- 
vant electric dipole transitions that are used to obtain the 
ion column densities from the data through Voigt-profilc 
fitting (see Paper I), are the same as used in the excitation 
modelling. 

3.1.1. Excitation modelling input flux spectrum 

The afterglow UV flux at the GRB facing side of the ab- 
sorbing cloud is obtained by converting the i?-band bright- 
ness as observed by the RAPTOR-T telescope (see the 
previous Sect. Kj) to the host-galaxy redshift at a partic- 
ular distance (ant parameter) from the GRB. This conver- 
sion includes a correction for both the Galactic extinction 



of Ar = 0.12 mag ( Schlegel et al 
ble extinction in the host galaxy. T 



1998 1, and any possi- 
le latter was found to 



be of type Small Magellanic Cloud (SMC), with an esti- 
mated U-band extinction (i n the host-galaxy rest frame) 
of Ay = 0.19 ± 0.05 mag ( |Kann et aL||2010[ ); we adopt 
this value for our main fits, and assume that the dust re- 
sponsible for this extinction is located within the absorbing 
cloud. The i?-band light curve is interpolated in log space 
to obtain the brightness at any given time to be used in 
the model calculations. To determine the flux at different 
frequencies, we initially adopt a value for the sp ectral slope 
of (3 = —0 .75 (the same as that adopted by Littlejohns 
et al.|2012 ). We note that /3 represents the intrinsic spectral 
jefore it is affected by the host-galaxy extinction 
The combination of this slope with an ex- 



slope, i.e. 
(if non-zero) 

tinction of Ay — 0.19 mag agrees well with the observed 
spectral slope: /3 = — 1, obtained from fitting the RAPTOR 
VRI data. Apart from this default slope-extinction setting, 
we also performed fits with zero host-galaxy extinction and 
the slope set to the observed value of (3 = —1. Since the 
afterglow flux in the X-ray regime is not relevant for exci- 
tation, we do not consider the X-ray flux. 

3.1.2. Excitation-only fit result 

The observed column density evolution of Fen, Fein, Sill 
and C II (both ground-state and excited levels) are fit with 
an excitation-only model, which includes the following fit 
parameters: 1) the GRB to cloud distance, i.e. the distance 
from the GRB to the front of the cloud, facing the GRB^J 
2) the linear cloud size, 3) the pre-burst Fen, Fein, Sin 
and C II column densities, and 4) the Doppler parameter 
describing the velocity distribution of the atoms. The re- 
sulting fit, shown in Fig. [2j describes the observed column 
densities quite poorly, with a large reduced chi-square value 
(x 2 = 21.3). One of the main reasons for the poor fit is that 
all observed levels of Fe II are decreasing with time, which 
cannot be accommodated with excitation alone. We note 
that - based on the atomic transition probabilities between 
the different levels - it is not possible for a large fraction of 
the pre-burst Fe n atoms to be excited to levels above those 
of the ground term ( 6 D). Moreover, transitions from these 
higher levels are not observed (see Paper I); e.g. see the Fen 
4 F 9 / 2 level upper limits (indicated by the blue triangles) in 



http: / /www. nist.gov/pml/data/asd.cfm 



Whenever we use the terms 'cloud distance', we refer to the 
distance from the GRB to the GRB-facing side of the absorbing 
cloud. 
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t rcst (minutes since GRB trigger) 




100 



t obs (minutes since GRB trigger) 



Fig. 2. Photo-excitation modelling of the observed (total) 
column densities as a function of time, as measured in 
Paper I (see their table 3), for Cn, Sill, Fen, and Fein. 
The open circles are detections, while the open triangles 
indicate upper or lower limits (3cr). The different colours 
denote the different ion levels: black for the ground state, 
and red-green-maroon-orange for the first four excited lev- 
els, while the Fe II 4 F 9 / 2 and Fe HI 7 S3 levels are indicated in 
blue. The model fit describes the observed column densities 
very poorly, with a reduced chi-square of xt = 21.3. 



Fig. [2j Another feature that is very difficult to explain with 
excitation alone, is the very large observed fraction of Fe in 
atoms in the excited 7 £3 level, of the order of 10% (see 
the blue level in the bottom panel of Fig. [2]) . This level is 
severely underestimated by the model fit, despite the best- 
fit cloud distance being lower than 50 pc. For these reasons, 
we can reject, with high confidence, the hypothesis that ex- 
citation alone is responsible for the observed column density 
evolution along the GRB 080310 sightline. 



3.2. Inclusion of photo-ionization 

Excitation of an ion in a particular ionization state does 
not change the total number of ions in that state. However, 
the observed Fen column densities all clearly decrease in 
time, including the ground state whose variability is gener- 
ally not detected. This suggests that Fell may be increas- 
ingly ionized (by the GRB afterglow) to higher ionization 
states such as Fein (see Paper I). This hypothesis is sup- 
ported by the detection of transitions of Fein, involving 
both the ground state and the 7 S3 excited level; this latter 
level has never been observed before along a GRB sightline. 
Moreover, in Paper I we find that the Fe ill 7 S3 level popula- 



tion is clearly increasing with time. Given this observational 
evidence for photo-ionization, and our finding above that 
photo-excitation alone cannot reproduce the column den- 
sity evolution observed, we have included photo-ionization 
in our model calculations. 



3.2.1. Modelling input flux spectrum: inclusion of X-rays 

With the inclusion of photo-ionization, we need to consider 
both the UV and X-ray afterglow radiation. X-ray photons 
photo- ionize species such as Fe II, Si II and Fe ill mainly via 
the ejection of inner-shell electrons. Since the X-ray flux 
for GRB 080310 is not a simple extrapolation of the opti- 
cal/UV flux with the optical spectral slope (see Littlejohns 



et al.||2012), we include the X-ray light curve as measured 



by the Swift X-ray Telescope (XRT). We retrieved the 0.3 
10 keV XRT afte rglow light curve in count rate from the 
Swift repository (Evans et al. 2009[ ) and separate it in 8 
different time intervals (0-141, 141-185, 185-269, 269-393, 
393-545, 545-615, 615-796, 796-7261 seconds after the trig- 
ger) in order to limit the possible X-ray spectral evolution 
in each single isolated light curve track. We then extract 
the spectrum from the repository for each time interval 
and we convert the count rate light curve to flux density 
accordingly. The monochromatic flux at 1.73 keV (logarith- 
mic average of the X-ray band) was calculated assuming 
the correspondent spectral slopes for each window. This X- 
ray light curve replaces the i?-band extrapolation in the 
regime above 0.3 keV (in the observer's frame, correspond- 
ing to 1.0 keV in the host galaxy rest frame). In the region 
0.3-10 keV we adopt the spectral slopes determined for 
the different time intervals, and beyond 1 keV we adopt a 
spect ral slope of f3 = —2 at all times (see Littlejohns et al. 



2012). The X-ray spectra and assumed spectral slopes are 
shown for two time intervals (185-269 and 615-796 seconds) 
in Fig. [3} 

We also performed fits with an alternative to the in- 
put flux spectrum described above. This alternative is 
motivated by modellin g of the GRB 080310 afterglow by 
|Littlejohns et al. (2012), which suggests that the early-time 
flux (up to about 1800 s in the observer's frame), in be- 
tween roughly 3 eV and 300 eV, is much higher than the 
j3 = —0.7 5 (or (3 = —1) extrapola tion from the optical (see 
figure 9 of Littlejohns et al.|2012 1 . We note that this regime 
has no observations that are able to constrain the proposed 
model. The Littlejohns model flux is approximated by in- 
terpolation of the RAPTOR optical and Swift X-ray light 
curves m between 3000 A (3.6 eV) and 300 eV (both in the 
observer's frame). Below and above this region the flux used 
is the same as the original input flux spectrum described 
above. In Fig. [3] we show the default input spectrum (solid 
line) and the Littlejohns alternative (dashed line), at two 
different epochs. In the modelling, the input spectrum is 
constructed by interpolation of the the RAPTOR i?-band 
and Swift XRT light curves for each new time step. 



3.2.2. Cross section of Fen ionization to different levels of 
Fe in 

Our progr am incorporates well-known astrophysical pro- 



log. 



Osterbrock & Ferland|[2006 ) , and we refer the 



reader to Appendix |A| for a detailed description of how the 
photons, with the input flux described above, that are trav- 
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^obs (A) 
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Energy obs (keV) 
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Fig. 3. The adopted input flux spectrum that is arriving at the GRB-facing side of the observed cloud - in observed 
flux units - depicted at two different epochs: t Q b s = 225 seconds on the left, and t bs = 672 seconds on the right. The 
solid line shows the default input spectrum with a spectral slope of (3 = —0.75 up to 0.3 keV, above which the X-ray 
flux (at 1.73 keV) and spectral slope is adopted. The RAPTOR VRI and Swift XRT observations corresponding to 
these epochs are overplotted with open squares. The long-dashed line shows the default input spectrum modified by an 
extinction of Ay — 0.19 mag; since this extinction is placed inside the observed cloud, it is not observable at the front of 
the cloud. The dotted line shows the input spectrum assuming a spectral slope of ft = — 1, combined with no extinction. 
The dashed line between the observed optical and X-ray regimes shows our approximation to the L ittlejohns model flux. 
The energy and flux limits of this figure are the same as those in figure 9 of Littlejohns et al. ( 2012 ) for easy comparison. 
The hashed regions show the flux decrease due to a foreground cloud with Hi column density log ./V(Hi) = 18.9 for the 
default input flux model (horizontal lines) and log N(H i) = 20.3 for the alternative Littlejohns input spectrum (vertical 
lines) - see Sect. Eland Table [3] for more details. The ionization edges of Hei at 24 eV/(l+z) and Hen at 54 eV/(l+z) 
can be spotted. We note that the spectral region in between the Lyman limit and the X-ray data is not constrained by 
imaging observations. 



elling through the cloud, are exciting and ionizing the ions 
therein. We stress that ionization is taken into account for 
all relevant ions, i.e. H i, He i, He II, Fe II, Fe in, Si II, C n and 
Cr II, and that we properly take into account the fraction 
of Fe II that will be ionized to Fe I II (rather than to higher 
ionization states), as calculated by Kaastra & Mewe (1993) 
for the different ion shells. Excitation is included for all 
ions except for hydrogen and helium. One very important 
non-standard aspect - the calculation of the cross section 
of Fe II ionization to different (excited) levels of Fe in - is 
discussed here. 

When Fen is ionized to Fein, the Fein ion will not 
necessarily be in its ground state, at least not immediately. 
We have calculated the photo-ionization cross section from 
Fen to specific levels of Fein, using two diff e rent c odes: 
the suite of programs developed by Cowan ( 1981 1 and 
the F lexible Atomic Code (FAC) developed by ]UuT 2003 



2004|. The Cowan code is a self consistent Hartree-Fock 
model with relativistic corrections. The FAC package is 
also a self-consistent program, which models the wave func- 
tions to self-consistency by including the electron screening. 
Relativistic effects are taken into account by means of the 
Dirac Coulomb Hamiltonian. 

The lowest configurations in Fen, 3d 7 and 3d 6 4s, 
strongly overlap, having 3d 6 ( 5 D)As 6 D g / 2 as its ground 
state. However the e D magnetic J-sublevels are just slightly 
higher in energy and are all populated (see Paper I). Since 
no absorption features from higher lying levels in Fe n have 



been observed (see table 3 of Paper I) , we have focused on 
ionization from the low lying 6 D levels. The character of 
the configuration that the ground state belongs to (3d 6 4s) 
results in a photo-ionization process that is both complex 
and interesting. There arc several channels to ionization 
from the ground configuration: (a) 3d 6 4s — > 3d 6 by ioniz- 
ing the outer 4s-electron to the p-continuum by means of a 
photon absorption, and (b) 3d 6 4s — > 3d 5 4s by ionizing the 
3d-electron to the p- or f-continuum by means of a photon 
absorption. 

The result of (a) is the population of the Fein 3d 6 5 D 
states and (b) will end up in the Fein 3d 5 ( 6 S)4s 7 S or 5 S 
states. In Fein there is quite an energy difference (3.7 eV) 
between 3d 6 5 D and 3d 5 ( 6 S)4s 7 S . The fact that absorption 
features arising from these states are observed and not from 
states in between, indicates that photo-ionization plays the 
dominant role rather than collisional excitation. 

In the approach using the Cowan Code the even con- 
figurations 3d 7 and 3d 6 As were applied and the odd con- 
tinuum states 3d 6 ep and 3e? 5 4s ep and ef. In case of FAC 
the Fen 3d 6 As and 3d 7 , and Fein 3d 5 4s and 3d 6 config- 
urations were introduced in the modelling. Comparison of 
the results from the Cowan and FAC programs shows a 
very good general agreement. Table [2] lists the calculated 
cross sections for the relevant levels of Fein and the cor- 
responding fraction of Fe n ionizations that populate that 
particular Fell! level. These numbers are used directly in 
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Table 2. Cross sections for ionization of ground-term Fell 
ions to different excited levels of Fein, as calculated with 
the FAC and Cowan codes. 



Fein level" 


cross section 


fraction of total 




(xlO" 19 cm" 2 ) 


% 


3d b »D A (1) 


3.50 


1.96 


5 D 3 (2) 


3.48 


1.95 


5 D 2 (3) 


3.47 


1.95 


6 -Di (4) 


3.46 


1.94 


5 A> (5) 


1.88 


1.05 


3d 5 (65')4s 7 S 3 (18) 


55.74 


31.26 


5 5 2 (26) 


1.69 


0.95 


3d b (4G)4s b G 6 (35) 


18.54 


10.40 


5 G 5 (36) 


19.03 


10.67 


5 G 4 (37) 


20.01 


11.22 


5 G 3 (38) 


22.49 


12.61 


5 G 2 (39) 


25.01 


14.03 



Notes. 



The Fell! level is indicated with the configuration, the term 



and subscript J value (and the level number — ordered in energy - starting 
from the ground level). 



our modelling programQWe find that only a small fraction 
(9%) will directly populate the Fein ground term, while 
31% of the new Fein ions will in fact populate the 7 S3 
level. The majority - 57% - will populate the levels of the 
3d 5 ( 4 G)4s 5 G term. However, these latter will quickly de- 
cay to the 3d 6 5 D ground term. Fig. [4] shows a partial en- 
ergy diagram of some relevant terms of Fe III. We note that 
many more terms exist, but are not depicted for clarity. For 
each relevant transition (indicated with the dotted line), 
we list the logarithm of the transition probability. E.g. for 
the strongest transition between the 3g? 5 4s 5 G and 3d 6 5 D 
terms, this is A = 10 +2 s _1 . The reciprocal of this num- 
ber provides the time in seconds in which the ions in the 
upper level would decay to the lower level in the absence of 
radiation. 



89084 



82001 



63425 - 



40999 - 



30088 



19404 



3d 6 




3d 5 4p 




3d 5 4s 



5 P° 2 



; +8.7 

.*' UV48 



• :> . ' +6.3 ; 

UV34: . 5 G 



1-8.6 



-5.3 : 




! 5 S 




-5.5 

7 S 



3 P2 




Fig. 4. Partial energy level (or Grotrian) diagram for the 
relevant lower terms of the three lowest configurations of 
Fein: 3d 6 , 3d 5 4s and 3d 5 4p (indicated at the top). The 
horizontal solid lines depict the energy levels, labelled with 
the term and J-value. Selected transitions are shown with 
dotted lines between levels. For each transition we list 
the logarithm of the transition probability (or Einstein A- 
coefficient, in [s -1 ]) of the strongest transition between 
the t erm s, adopting the Raassen & Uylings values (see 
Sect.O). 



3.2.3. Comparison with ionization calculations in the 
literature 

As a consistency check, we compared the amount of ion- 
ization computed in our program with calculations in the 
literature. Since our ionization model does not take into 
account recombination^] while most calculations in the lit- 
erature do, the comparison options are limited to GRB ion- 
ization studies. E x amples of these are the s tudies of Perna 
fe Lazzatil ( |2002| ), |Perna et ah] ( |2003[ ) and |Draine fe Hao 
( 2002J, in which not only the ionization induced by the 
GRB is calculated, but also the accompanying destruction 
of dust and dissociation of H 2 . Since our program does 
not include dust destruction, this makes a comparison with 
these calculations difficult. However, we were able to com- 
pare our progra m with the H, He , and N photo- ionization 
calculations by Prochaska et al. (2008), and find consis- 
tent results. We use their Eq. 8 torTne GRB 050730 after- 



7 We note that before we had calculated these numbers, we 
included the fraction of Fen ionizations that populate the 7 S3 
level of Fe ill as a free parameter in our model fit, with a resulting 
best-fit value of 30-35%. 

8 This is not required, because the relevant time scale of our 
calculations - hours to days after the GRB - is negligible com- 
pared to the recombination time scale at typical ISM densities. 



glow luminosity over the same time span i ODS = 10-1000 s 
and adopt their set-up with tt-h = 10 cm~ 3 , a nitrogen- 
to-hydrogen abundance ratio of 10 -6 (roughly 0.01 solar 
metallicity) and assume that before the GRB all the ions 
are in the singly ionized state. We then switch on the 
GRB 050730 afterglow and follow the progressive ionization 
of N 11 to higher ionization states, and find that after 1000 s, 
the Nv column density remaining is log 7V(Nv) = 13.8, 
compared to their log A r (Nv) = 14. Also the ionization 
structure at £ ODS = 1000 s computed by our program is 
very similar to that depicted in their figure 3. 



3.2.4. x 2 Minimization and fit parameters 

The model column densities computed by our program, as 
detailed in Appendix [A] are fit to the observed GRB 080310 
column densities at their respective epochs. We use the 
Fortran 90 ve rsion of the MINPACK lmdif \ 2 minimiza- 
tion routine (More et al. 1984 1980), which is based on 



the Levenberg-Marquardt method. We have made this pro- 
gram parallel with OpcnMP, so it can be run faster on a 
shared-memory computer cluster. The formal errors of the 
fit parameters are estimated by computing the co-variance 
matrix, and taking the square root of the diagonal elements. 
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Fig. 5. Photo-excitation and -ionization modelling of the 
observed (tot al) column densities as a function of time, as 
measured by De Cia et al. (2012) (see table 3 of Paper I), 
for - from top to bottom panels - Hi, Cn, Sin, Cm, Fen, 
and Fein. The solid and dotted lines correspond to the 
best-fitting model assuming the default and the Littlejohns 
input flux, respectively. The different colours of the symbols 
and lines have the same meaning as in Fig. [2] Although 
overall the model provides a reasonable description of the 
column density evolution, the observed excited levels of Fen 
at epoch II are significantly underestimated. See the text 
and Table [3] for more details. 



The fit parameters are the same as those used in the 
excitation-only case described at the end of Sect. |3.1[ the 
GRB to cloud distance (i.e. the distance from the GRB 
to the GRB-facing side of the cloud), the cloud size, the 
Doppler parameter 6, and a pre-burst column density for 
each ion included in the fit. We again initially fix the spec- 
tral slope to j3 = —0.75 (see Littlejohns et al. 2012 ), in com- 
bination with a h ost-galaxy extinction of Ay = 0.19 mag 
( Kann et al.||2010 1 , but also experiment with the combina- 
tion fj = —1.0 and zero extinction. The ions included are: 
Fen, Fein, Sin, Cn and Cm (apart from Hi, Hei and 
Hen, see below), i.e. all low-ionization species with a total 
column density measurement at one or more epochs as re- 
ported in table 3 of Paper I. For all of these, except Cr II, 
we include excitation. 

As discussed in Paper I, the velocity profiles of Fen and 
Fein are markedly different from those of Sin and C II. The 



60 
O 




t obs (minutes since GRB trigger) 



former are dominated by component "b" at —20 km s 



100 



Fig. 6. Same as Fig. |5j but including a foreground cloud 
situated in between the GRB and the observed absorber. 
The solid and dotted lines correspond to the best-fitting 
model assuming the default and Littlejohns input flux, re- 
spectively. The foreground absorber is mostly ionized by 
the time of the 2nd epoch UVES spectrum and can there- 
fore escape a clear detection. See the text and Table [3] for 
more details. 



from the systematic velocity and with a Doppler broadening 
parameter value of 6 = 13 km s _1 . However, a considerable 
column density is contained in the other components "a", 
"c" and "d" as well, leading to an overall broad velocity 
structure for Fen and Fein. In contrast, the vast major- 
ity (although not all) of the Si 11 and C II column densities 
is located in the narrow "c" component at the systematic 
velocity, with 6 = 7 km s _1 . For this reason, we split the 
velocity broadening fit parameter into two: one for Fe 11 and 
Fein: fcpeii.FeiiD an d one f° r Sin and Cn: 6sin.Cn- 

In Paper I, we also constrained log A^(H 1) to 18.7 at two 
different epochs. This H 1 column density is an important 
quantity, as it can - if it is sufficiently large effectively 
shield the low-ionization species (such as Fen and Sin) 
from the ionizing photons. Besides Hi, we also include Hei 
and He 11, which are also important for shielding, albeit at 
higher photon energies (starting from 24 eV) . These helium 
ions do not require additional fit parameters, as we fix the 
He I colu mn density at the sola r abundance value (i.e. 8.5% 
of N(H 1), [Asplund et aL]|2009| ) and set the pre-burst Hen 
column density to zero. We note that the inclusion of a 
significantly larger amount of He 1 - a possibility if there is 
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Table 3. Column-density evolution modelling fit results. 



Incl. foreground cloud? no 


yes 


yes 


Incl. Littlejohns flux? no 

2 1 i r r i \ „ - i 


no 


yes 



Cloud distance (pc) 
Cloud size (pc) 
&Feii,Fem (km s _1 ) 
&Sin,Cn (km s" 1 ) 
log JV(Hl) [cm" 2 ] 
log AT (Fell) [cm -2 ] 
log JV(Fem) [cm -2 ] 
log JV(Sill) [cm" 2 ] 
log iV(Cll) [cm" 2 ] 
log JV(Crn) [cm" 2 ] 
FC d distance 6 (pc) 
FC d size (pc) 
FC d log AT (Hi) [cm" 



363 ± 86 
0± 149 
50 c 

2.f ±0.7 
18.75tg : » 



14.. 



1+0.03 
J -0.03 

14-56±g;?g 
13.70l° ° 6 7 
14.28l° °° 
13.52±?," 



235 ± 97 
126 ± 141 
50 c 

1.5 ±0.6 
18.64l° ;l 2 7 
14.831°;°! 
14.60+ 008 
13.74 
14 28+°'° 6 

13 49 +0.20 

50 c 

39 ± 715 
18.9±°;g 



0.10 
0.07 
0.0K 



-0.37 



260 ± 107 
55 ± 146 
38 ± 13 
1.4 ±0.5 
18.60 
14.93 
14.70 
13 

14.30 
13.60 
12 ±8 
9± 19 
20.30±g-jS 



+0.11 

0.14 
+0.07 

0.08 
+ 0.10 

0.13 
+0.07 
-0.09 
+0.07 

0.08 
+0.13 

0.18 



Notes. 



(a) 



The spectral slope and host-galaxy extinction were fixed to 



the values: /9 — —0.75 and Av — 0.19 mag, respectively. ^ This is the 
distance from the GRB to the GRB-facing side of the cloud. ^ Maximum 
allowed fit value reached. ^ Foreground Cloud. 



a large column density of pre-burst ionized hydrogen - does 
not affect our results. If the Hei abundance is included in 
the fit as a free parameter, the best-fit value is consistent 
with the adopted value: [Hel/Hi] = (12 ± 17)%. 



4. Results 

The resulting fit to the total column densities is shown 
in Fig. [5j The solid (dotted) lines correspond to the best- 
fitting model assuming the default (Littlejohns) input flux 
discussed in Sect. |3.2| The goodness-of-fit and best-fit pa- 
rameter values are listed in the first column of Table [3] 
The model fit in which the Littlejohns input flux is adopted 
is very poor, with xl — 8.34, and we therefore discard it 
without listing the unreliable best-fit parameter values in 
Table [3j The quality of the model fit that assumes the de- 
fault input flux is reasonable, with xl — 2.51. This rather 
high value for the reduced chi-square seems to be mainly 
caused by the model underpredicting the observed pop- 
ulation of the Fen excited levels. Assuming a negligible 
host-galaxy extinction (Ay = mag), combined with the 
observed spectral slope of /3 = — 1, leads to a slightly im- 
proved fit with a chi-square value of xl — 2.36, but with 
resulting best-fit values consistent within the errors of the 
default fit (with (3 = -0.75 and Ay = 0.19 mag). 

Table [3] shows that the best- fit Si 11 and C 11 Doppler 
broadening parameter is very low: fosin.Cn = 2.1 ± 0.7 km 
s _1 . As we discussed in the previous section, the observed 
6-parameter value is low as well: &sin,Cn — 7 km s . To 
investigate this modest discrepancy further, we also ran a 
model in which only Hi, He I, Sin and C11 are included, 
i.e. without Fell, Fein and Cm, and with the ^-parameter 
fixed to the observed value of b = 7 km s _1 . Although the 
resulting distance to the GRB-facing side of the cloud is 
very small, less than 50 pc, the cloud size becomes more 
than a kiloparsec, i.e. the average distance is quite large. 
Indeed, forcing the cloud to be very compact, with a cloud 
size fixed at 1 pc, the best-fit distance, both without and 
with the Littlejohns input flux, is roughly 600 pc. These re- 
sults indicate that the majority of Si 11 and C II ions might 



be at a different location (further away from the GRB) or 
spread out over a larger region than the bulk of the Fell 
and Fe ill ions. This is supported by the very different veloc- 
ity profiles that these ions display (see Paper I). However, 
for other GRB sightlines for which a cloud distance has 
been determined independently for Fe II and Si 11 excitation 



(e.g. D'Elia et al. 2010 2011), the best-fit distances are 



consistent. This suggests that - as one would expect - the 
Fen and Sin atoms are probably located at comparable 
distances from the GRB. 

Since the Fe 11 excited levels are under-predicted by the 
model, we attempted to place an additional cloud along 
the line of sight, in between the GRB and the observed ab- 
sorber or cloud. If the additional cloud would be sufficiently 
close to the burst, it would become completely ionized dur- 
ing the first few tens of minutes (in the observer's frame) 
of the arrival of the GRB radiation, and would not reveal 
. itself in the observations. But at the same time it would 
partially shield the observed cloud from ionizing radiation 
released during the first minutes after the GRB, allowing 
the observed cloud to be closer to the burst, increasing the 
amount of excitation. This scenario could only work in case 
the two clouds have a velocity offset (already 10-20 km s _1 
is sufficient, which is not unlikely), to avoid the observed 
cloud to be in the absorption-line shadow of the foreground 
cloud. The fit with such an additional cloud is shown in 
Fig. [6j again with the solid (dotted) curves corresponding 
to the model fit adopting the default (Littlejohns) input 
flux, and the best-fit parameter values are listed in the 2nd 
and 3rd column of Table OH The addition of such a fore- 
ground cloud results in a lower value for the chi-square: 
xl = 1-50 (assuming the default input flux), but at the 
expense of three additional fit parameters: the distance, 
size and column density of the foreground cloud (FC, see 
Table [3]) . An F-test suggests that the fit improvement in- 
troduced by the foreground cloud is significant, providing 
F = (x 2 -x|c)/(^fc) = 4>8 ( w here v is the number of de- 

Xpc/fFC V 

grees of freedom) and a null probability of P < 0.005; i.e. 
there is less than 0.5% chance that such an improvement is 
random. 

We also ran model fits with both a foreground cloud 
and adopting th e alt ernative Littlejohns input flux. As de- 
this input flux is much higher (up 



scribed in Sect. 3.2 



to a factor of 10) than the default input flux i n between 
0.3 eV and 300 eV (see |Littlejohns et al ||2012| ), i.e. lead- 
ing to much more ionizing radiation. As already mentioned 
above, a model fit with the Littlejohns input flux without 
an additional cloud describes the observed column density 
evolution very poorly. However, an additional cloud with a 
neutral hydrogen column density almost that of a damped 
Lya (DLA) system at 10-20 pc from the GRB, is capable 
of absorbing most of the extra ionizing radiation, leading 
to a very reasonable fit (with xl = 1-71). The best-fit pa- 
rameter values for this model are listed in the 3rd column 
of Table [3j and the resulting column-density evolution is 
shown with a dotted line in Fig. [6j 



5. Discussion 

The evolution of the Fen and Fein column densities 
observ ed a t the GRB 080310 redshift (see Paper I, and 
Sects. 



3.2 and |4j), combined with our modelling, clearly 
shows that ionization of Fe 11 is taking place. A very strong 
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argument in favour of ionization and a vital ingredient for 
the model ling, is that according to our calculations (see 
Sect 3.2.21 a large fraction (31%) of Fen ionizations will 
- initially - populate the Fein 7 S3 level. Without taking 
this effect into account, we found it impossible to explain 
the large fraction (~10%) of Fein that is observed to be 
in this particular level (see Paper I). This channel for pro- 
ducing a significant Fein 7 S3 level population may be rel- 
evant for other objects in which absorption lines from this 
level - the UV34 triplet - are also observed, such as broad 
absorption line (BAL) quasars and r\ Carinae. As it takes 
about 1000 seconds for the Fein 7 53 level population to 
decay spontaneously down to the ground term, the Fell 
ionization rate needs to be significant at this time scale 
for this process to be relevant. In BAL quasars, the UV48 
triplet, at 2062, 2068, and 2079 A, is sometimes detected. 
The lower energy level from which the UV48 triplet arises, 
5 ^2 (see Fig. [4K, receives only a small fraction of the Fen 
ions that are ionized to Fein (1%, see Table [2]), and so 
these lines are expected to be much weaker than the UV34 
triplet. We checked for the presence of these UV48 absorp- 
tion lines in the UVES spectra of GRB 080310, and indeed 
do not detect them. In the sample of unusual BAL quasars 
of |Hall et al] ( |2002[ ), the detection of the UV34 triplet is 
much more common than UV48, which would be expected 
in case ionization of Fe 11 is the dominant mode of popu- 
lating the UV34 lower level. However, in case the UV48 
absorption is stronger than that of UV34, such as in SDSS 



2215-0045 (Hall et al. 2002 Vivek et al. 2012), the above 



Fe in excited- level population scenario, which works well for 
GRB 080310, does not provide a viable explanation. 

Time- variation of H 1 and metal-column densities due to 
the ongoing ionization by the GRB and afte rglow radiation 
has been predicted (e.g. Perna fe Loeb | 1 99 8]), but ha s never 

"u — : 1_- \ 1 1 — r / ^ ^^T^f77~ „l Horn 1 V 



been convincingly detected before (see Tho ne et a,l.|[20lT b. 
This applies not only to neutral-medium ions such as H 1 
and Feu, but also to high-ionization spec ies as C iv and N V 



(Prochaska et al. 2008 Fox et al. 2008). The reason that 



ongoing photo-ionization is observed for GRB 080310 is not 
that the observed neutral material along the GRB 080310 
sightline is much closer to the GRB than in other cases. 
We find a distance range of 200-400 pc (depending on 
the adopted input flux and the inclusion or not of a fore- 
ground cloud, see Sect. [4] and below), while other GRBs for 
which only excita tion was detected h ave distance estimates 
as low as 50 pc ( D'Elia et al.||2011[ ). In Table [I] we have 
collected the GRB-cloud distance estimates from the liter- 
ature, allowing for a direct comparison with the distance 
estimate for GRB 080310. We note that, in the absence 
of a foreground cloud, a lower limit of about 100 pc can 
be placed on the GRB-absorber distance, by just consid- 
ering the non- variation in the Hi column density between 
21 and 50 minutes post-burst. If the absorber would have 
been much closer, we would have detected a significant H 1 
column-density change. 

We investigated whether the very low Hi column den- 
sity or super-solar iron abundance along the GRB 080310 
sightline is the reason for the unique detection of ongoing 
ionization. We did so by running models with most pa- 
rameters fixed to the best-fitting model (with xl = 1-5 m 
Table[3]), but varying the H 1 column density and iron abun- 
dance, to see how these affect the number of ions detected 
in the Fe ill 7 S3 level. As we have shown above, a significant 
population of this excited level is a clear sign of ongoing ion- 



ization of Fen. Table [5] shows the expected peak Fein 7 S3 
column density as a function of different H 1 column densi- 
ties (rows) and iron abundances (columns). In the column 
with fixed -/VfcII.FcIH, we fixed the pre-burst Fen and Fein 
column densities at their best-fit values (middle column) of 
Table [3] for the four different Hi column densities, while in 
the last two columns, all the iron was assumed to be in the 
singly ionized state before the GRB exploded; we note that 
this latter assumption is only valid at higher H 1 column 
densities (log iV(Hi) > 20). 

Table 5. Maximum Fein 7 S3 column density reached as a 
function of the assumed H 1 column and iron abundance in 
the GRB 080310 absorber. 



log iV(Hi) 


fixed A^ eII , Pe „i 


[Fe/H]= 


-1.0" [Fe/H]= +0.2"< c 


18.6 


13.55 


11.83 


13.02 


19.6 


13.24 


12.52 


13.71 


20.6 


12.69 


12.95 


14.14 


21.6 


12.05 


13.20 


14.37 



Notes. In this column the Fell and Fein column densities were fixed 
at the best-fit values of Table [3| i.e. for higher H I column densities the iron 
abundance effectively decrease with respect to the log iV(Hl) — 18.6 case. 
^ For these runs at fixed iron abundance, we have assumed that before 
the onset of the GRB, [Fc/H] — [Fell/Hi], i.e. all the iron ions are in the 
singly ionized state and all the hydrogen is neutral, which is only a good 
approximation at higher Hi column densities. ^ The iron abundance for 



GRB 080310 was determined to be [Fc/H]: 
corrections (sec Paper I). 



=+0.2 when including ionization 



Considering log iV(Felll 7 S 3 ) = 13 to be the approxi- 
mate lower limit for a clear detection of this level in the 
UVES spectra, Table [5] shows that increasing the H I col- 
umn density by a factor of about 100 or more, while fixing 
the Fen and Fein column densities at the best-fit values 
of Table |3j would have resulted in a non-detection of the 
Fein excited level. This is due to the increased shielding 
of the low- ionization metals from the ionizing radiation by 
the Hi and Hei atoms. However, when fixing the abun- 
dance at the observed value for iron along the GRB 080310 
sightline: [Fe/H]=+0.2 (see Paper I), the Fein excited level 
is detected at any Hi column density. At a more typi- 
cal iron abundance for GRB sightlines, [Fe/H]= —1.0, the 
Fein UV34 triplet is detectable only at the higher Hi col- 
umn density end. We note that a column density of log 
A r (Hi) = 21.6 at 0.1 Z Q implies a considerable Fell col- 
umn: log AT (Fen) = 16.1, and this increases by at least a 
factor of 10 when assuming [Fc/H]= +0.2. At such large 
Fen columns, dust is not unlikely to be present. The pres- 
ence of dust would complicate the UV34 triplet detection 
at high Hi columns. Dust obscuration would not only de- 
crease the amount of ionization taking place, but would 
also make it more difficult to detect a bright afterglow - 
required to secure high-quality spectra - in the first place. 
Therefore, the reason for the unique detection of the Fe in 
UV34 triplet in the GRB 080310 spectra appears to be a 
combination of the super-solar iron abundance and the low 
Hi column along this sightline. This ensures the presence of 
a sufficient amount of iron, while at the same time avoiding 
too much H I and He I shielding and dust obscuration. 

If the detection of Fen ionization is indeed due to a 
combination of the super-solar iron abundance and the low 
H I column density along the GRB 080310 sightline in the 
host, then the (non-)detection of the Fe III UV34 triplet can 
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Table 4. GRB absorber distances, as of April 2012. 



GRB 


Instrument 


z 


Distance 

(PC) 


Size 

(PC) 


pa 


log AT (Hi) 


[X/HJ 


X 


Ref. 


020813 


LRIS+UVES 


1.25 


50-100 












1, 2 


050730 


UVES 


3.97 


124 ± 20 6 


147Ij4 


E 


22.10 ± 0.10 


—2.18 ± 0.11 


S 


3, 4 


Uollll 


TJT"D T^G 

rilruib 


1 K K 

1.00 


a tew nunurea 




j? 
hj 








0, O 


060418 


UVES 


1.49 


480 ± 56 




E 


> 21.0 


< -0.5 


Zn 


7, 4 


080310 


UVES 


2.43 


200-400 


0-200 


E+I 


18.70 ±0.10 


-1.2±0.2 


Si 


8, 9 


080319B 


UVES 


0.94 


560-1700 




E 








10, 4 


080330 


UVES 


1.51 


7 q+il 




E 








11, 4 


081008 


UVES+FORS 


1.97 


52±6 C 




E 


21.11 ± 0.10 


-0.87 ±0.10 


Si 


12 


090426 


LRIS+FORS 


2.61 


£80 d 




I 


18-7±g-i 






13 


090926 


X-shooter 


2.11 


677 ± 42 e 




E 


21.60 ±0.07 


-1.85 ±0.10 


s 


14, 4 



Notes. The distances derived from a photo-excitation/photo-ionization model of the column density variability, based on h igh-resolution spectrosc opy, 
highlighted in bold, arc considered to be more reliable. The excitation distances have been corrected for yj 4n according to | Vreeswijk et al. 
^ P rocess modelled: photo-excitation (E) o r photo-ionization (I). ^ A former analysis of the Magellan Clay/MIKE echelle spectrum " |ch en et al. 
2005 1 suggested a cloud distance d < 100 pc l Prochaska et al.|2006 l. (c> Component I (Component II lies at 200+™ pc). (d) We consider the 090426 
Lya variation detection to be marginal, and conservatively list this distance estimate as a lower limit. ^ Main component (the second lies at ~ 5 kpc). 

Reference s. ( 1) |Dessauges-Zavadsky et al. 1)20061; (2)|Savaglio fc Fall l (|2004| >; (3) |Ledoux et al.|||2009[; (4 ) Vreeswijk et al.| j20lT|> ; (5) |Penprase | 
let al.|(|2006h (6) |Prochaska et al.|liM)fl[; ( 7) | Vreeswijk et al.|0uO7t; (ti ) Thi s work; (9) Paper I: ILe (Jla et ah] ; (10) llVUl'a. et al.| (2M{>a[; (11) 

|lVll!l!a et af] ( |20(fob[ ; (12) |U'll!l!a et al.| ( [201 1| ; (Uj) |fhone et al.| ( |201l[ ; (14) |D'Elia et al.|p010[ 



be used to put constraints on the H I column density along 
a GRB sightline with an Fe II detection in case it cannot be 
inferred from the spectrum. The latter is the case at z 2, 
when Lya is not redshifted enough to be included in the 
optical wavelength range of spectrographs on ground-based 
telescopes. The strength of the Fein UV34 triplet, however, 
depends on various quantities besides the iron abundance 
and Hi column, such as the GRB-absorber distance, the 
afterglow peak luminosity and brightness evolution, and the 
time at which the spectra are taken. It is therefore difficult 
to provide a simple scaling relation between the H I column, 
iron abundance and UV34 triplet strength. 

But for GRBs for which most of the above quantities can 
be constrained through absorption-line photo-excitation 
modelling, it is possible to determine a lower limit on the H I 
column density from the Fe in UV34 triplet non-detection. 
As our team has already p erformed such model l ing on 
GRB 060418 (at z = 1.490, |Vreeswijk et al.||"2007| |2011| , 
we can readily determine this limit on H I for this sightline. 
The UV34 triplet is not detected in the GRB 060418 UVES 
spectra, with a 3a upper limit on the rest-frame equiva- 
lent width (column density) of 0.03 A (log jV(Hl) = 12.6). 
Modelling the excitation and ionization with our code, in 
which we vary the H I column density, we find that this 
UV34 detection limit corresponds to an H I column density 
limit of log N(H i) > 21.0. Using the total zinc column den- 
sity measured for this sightline (log TV(Znii) = 13.09±0.01, 
Vreeswijk et al. 2007), and assuming that most of the zinc 
is in the singly ionized state, the Hi column density lower 
limit derived above implies an upper limit on the metallic- 
ity of [Zn/H]< —0.5. Determining these Hi column-density 
and corresponding metallicity limits for the entire sample 
of Table [4] requires photo-excitation and -ionization mod- 
elling for each sightline separately, which is out of the scope 
of the current paper. 

Our simplest model, in which the GRB afterglow is ion- 
izing and exciting a cloud at a distance of about 360 pc, 
does not provide a satisfactory description of the observa- 
tions. As can be seen in Fig. [5j the model underestimates 
the ground-term fine-structure level population. One poten- 



tial reason for this lack of Fe II excitation - or abundance 
of ionization - may be that additional neutral material is 
present in between the GRB and the absorber responsible 
for the absorption features observed in the spectra. This 
additional absorber needs to be ionized by the time that 
the first couple of spectra were taken, as otherwise it would 
reveal itself in the observed spectra. Placing such an addi- 
tional cloud closer to the GRB, with log 7V(Hi) ~ 19 at a 
distance of tens of parsecs, improves the model fit signifi- 
cantly, as shown by the solid curves in Fig. [6] Also in the 
case where the Littlejohns input flux is adopted (depicted 
by the dotted curves in Fig.|6|, the model with a foreground 
cloud provides a very reasonable description of the observed 
column-density evolution of the different ionic species. In 
this case, the foreground cloud is required to have a higher 
neutral hydrogen column density (log iV(Hi) = 20.3) and 
to be closer to the GRB (12 pc), in order to be able to 
absorb the additional ionizing photons in the Littlejohns 
input flux. The similar chi-squares for the additional-cloud 
model using the default and Littlejohns input fluxes does 
not allow us to favour one input flux over the other; how- 
ever, in the model without an additional cloud, the default 
input flux is clearly favoured. 

We tested if a log 2V(Hl) = 20.3 cloud at 12 pc, with 
an assumed metallicity of one tenth of solar and using the 
Littlejoh ns input flux, would imply an observab le N V varia- 
tion (see |Prochaska et al.|2008{|Fox et al |2008| ) in our spec- 
tra. In Paper I we report a constant Nv column density: 
log 7V(Nv) = 14.10 ± 0.04 and log 7V(Nv) = 14.05 ± 0.02 
at epochs II and IV, respectively. In this test, we adopt a 
metallicity of one-tenth of the solar abundance, and we as- 
sume that all the nitrogen is singly ionized before the burst. 
We find that the N n ions are very quickly ionized to higher 
ionization states: at 6 minutes after the burst (observer's 
frame) the Nv column density in the foreground cloud is 
already below log 7V(Nv) = 13, and by the time of the 1st 
epoch spectrum (13 minutes after the burst), practically 
all the nitrogen has been ionized to states higher than N v. 
Also, if the foreground cloud is indeed ionized within about 
ten minutes of the arrival of the first GRB photons, it is 
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very difficult to infer its presence in sightlines where only 
ongoing excitation is observed. 

Although the introduction of an additional cloud is a 
rather ad-hoc solution for improving the model fit of Fig. [5]- 
mitigating the ionization with respect to the amount of ex- 
citation - the existence of an additional cloud in the vicinity 
of the burst is not unexpected, as GRBs are tho ught to oc 



cur in gas- rich massive-star forming regions (e.g. Prochaska 



et al.|2007 ) . We note that the presence of such an additional 
cloud is consistent with the host-galaxy iV(H)-equivalent 
X-ray absorption as inferred from the Swift XRT data 
(logiV(H) = 21.7 ± 0.05 and log N(R) < 21 - assuming so- 
lar metallicity - for the time-averaged averaged Windowed 
Timing and Photon Counting modes, respectively, |Evans| 



et al. 20091. Thus, although the presence of a foreground 
cloud is plausible, we cannot exclude a different origin for 
the underestimate of the Fen excitation (or overestimate 
of ionization) in our default model fit. 



6. Conclusions 

We modelled the variability of the ionic column densities of 
various species (including Hi, He I, Hen, Fen, Fein, Sin, 
C II and Cr n) in the circumburst medium of GRB 08 0310 
(reported in a companion paper by: De Cia et al. 2012 ) with 
a photo-excitation and -ionization radiative transfer code. 
The rest-frame near-infrared to X-ray spectrum of the after- 
glow radiation and its time evolution, an important input 
parameter in the modelling, is inferred by combining the 
RAPTOR-T VRI light curves - also presented in this pa- 
per - and the X-ray light curve as observed by Swift. We 
find that excitation alone, which has been successfully ap- 
plied to other GRBs, is not able to explain the GRB 080310 
observations; ionization is clearly required. The strongest 
evidence for ionization is presented by the clear detection 
of the UV34 triplet of Fein from the lower level 7 S3. The 
large fraction of Fe in ions measured to be in this level, 10%, 
can only be explained through ionization of Fe II; we calcu- 
late that 31% of all Fen ions that end up as Fein will first 
populate this 7 S 3 level. This is the first conclusive evidence 
for the detection of time-variable photo-ionization induced 
by a GRB afterglow. 

Despite this evidence for photo-ionization, the distance 
between the GRB and the absorbing medium that we infer 
(200-400 pc) is very similar to that in other GRB sight- 
lines for which such a distance estimate was possible. We 
find that the main reason for detecting time-variable ion- 
ization in this GRB and not in others is the super-solar iron 
abundance ([Fc/H]=+0.2) in combination with the low Hi 
column density (log iV(Hi) = 18.7 ± 0.1) along this sight- 
line. 

The combined photo-excitation and -ionization mod- 
elling provides tentative evidence for the presence of an 
additional absorbing cloud, with log iV(Hi) ~ 19-20, at 
a distance of 10-50 pc from the GRB, even though this 
cloud is almost completely ionized by the afterglow within 
a few tens of minutes (in the observer's frame) of the arrival 
of the GRB radiation. Future time-resolved high-resolution 
spectroscopic observations of low-H I GRB sightlines could 
provide additional constraints on the existence of pre-burst 
neutral gas in the GRB vicinity. 
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Appendix A: Details of the time-dependent 
photo-excitation and -ionization calculations 

This appendix describes in detail our time-dependent 
photo-excitation and -ionization calculations of the neutral 
medium nearby the GRB, along the line-of-sight. Although 
this includes well-known astrophysical processes, it allows 
for a transparent comparison with similar future studies. 

Our program is rather basic when compared to a photo- 
ionization code such as CLOUDY. It does not include 
recombination, which is a reasonable assumption due to 
the very short time scale (of the order of hours to a 
day) that the GRB afterglow is bright and that high- 
resolution spectra can be secured. With an approximate 
rate of 10~ 13 cm 3 s _1 , the recombination time scale at typ- 
ical ISM densities is orders of magnitude larger. Second, 
our ID calculations are performed along the line-of-sight 
only and we do not take into account afterglow photons 
that have scattered off particles elsewhere in the absorb- 
ing medium and into the sightline. And third, we consider 
the source of afterglow photons to be very small compared 
to the distance from the source to the absorbing medium. 
However, the important asset of our code is its ability to 
take in a time-variable input source, and calculate the re- 
sulting column density evolution, for both the ground state 
as well as excited levels, as a function of time. 

Photo-excitation and -ionization are treated separately 
but self-consistently in our model. Photo-excitation in- 



volves bound-bound transitions in the ion, while photo- 
ionization involves bound-free transitions; hence the excit- 
ing photons have quasi-discrete wavelengths, while the ion- 
izing photons constitute a continuum with photon energies 
larger than the ionization threshold. Therefore, for the ex- 
citation we use a flux array that is calculated at specific 
wavelengths only - the central wavelengths of the relevant 
transitions - while for ionization we use a continuum flux 
array. This continuum array starts at the lowest ionization 
threshold of the ions used in the calculation (typically Hi 
at 13.6 eV) up to 100 keV, with logarithmic increments to 
have most of the wavelength resolution around the lower 
photon energies. Ideally, we would use a single continuum 
flux array, and calculate the imprint of both excitation and 
ionization as the photons move into the absorbing cloud, 
but this would be very CPU intensive as such an array 
would require both a large wavelength range, from the far- 
infrared to hard X-rays, as well as a sufficient resolution 
over this range. The two flux arrays overlap, but over a 
fairly short wavelength range. E.g. excitation of Fen ions 
is due to photons at particular wavelengths up to the ion- 
ization potential of Fen, i.e. 16.2 eV, and therefore they 
are in the range of the ionizing flux array. In this overlap- 
ping region, we take into account the flux decrease due to 
ionization on the exciting flux array. 

The absorbing cloud is divided into a number of plane- 
parallel layers. The initial (pre-burst) column densities of 
the cloud are spread out evenly between the layers such 
that each layer is composed of sublayers of pure Hi, Hei, 
Hen, and any additional ion such as Fell, Fein, and Sin. 
The number of layers is determined by our requirement 
that each layer is optically thin (r < 0.05) both at the Hi 
ionization threshold, as well as at the central wavelength of 
any excitation transition. To mimic a continuous excitation 
and ionization proces, we also require not more than 5% of 
the ions in any sublayer to be excited or ionized in any 
iteration, or time step. This - together wi th t he amount of 
flux arriving at the cloud front (see Sect. |3.1| and below) - 
determines the time step, At, of the calculation. 

The flux, F v (in erg cm~ 2 s _1 Hz" 1 ), is known at any 
given time and distance from the GRB. For the spectral 
range from the infrared to 300 eV, we use the RAPTOR- 
T i?-band light curve and adopt a spectral slope of (3 = 
—0.75 or f3 = —1, depending on the assumed amount of 
host-galaxy extinction (see below). Above 300 eV we use 
the Swift X-ray light curve, as described in Sect. |3.2| The 
observed fluxes are converted to the host-galaxy rest frame. 
For each time step, starting at 25 s in the rest frame of 
GRB 080310, up to the epoch of the last UVES spectrum, 
the corresponding flux is propagated through the cloud. 
It decreases as F n i ayor = i^iaycr- 1 e~ r - with r being the 
optical depth - as in every (sub-)layer photons are exciting 
and ionizing ions. The propagating flux is also affected by 
the dust extinction in the absorption med ium, which is fixed 
to Ay — 0.19 mag (Kann et al. 2010) when adopting a 
spectral slope of /? = —0.75, or to zero m combination with 
P = — 1. The extinction at wavelengths other than at the 
centre ofthe_K-band filter is calculated using the analytic 
fits of Pei (1992), assuming an SMC-type extinction. 



For the excitation optical depth, we use: 
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where H (a, u) is the Voigt function, a and u being func- 
tions of the frequency v and the Doppler width b. e and m e 
are the charge and mass of the electron, respectively, / is 
the oscillator strength, A the central rest-frame wavelength 
of the transition, N the column density, g the statistical 
weight of the level, and subscripts 1 and u indicate the lower 
and upper level, respectively, of the transition in question. 
The flux arriving at a particular layer is determined at the 
line centre of the transition. As long as the absorption line is 
still in the optically thin regime, this is equal to integrating 
the flux over the line profile. However, the line-centre flux 
starts to diverge from the profile-averaged flux at larger 
optical depths. The ratio of these fluxes depends on the 
optical depth only (and not, as one might expect, on the b- 
parameter of the absorbing medium) , and therefore we can 
easily convert the flux at line centre to a profile-averaged 
flux at any optical depth through tabulation of this ratio 
at a range of optical depths. 

When considering a cloud with a non-zero size, the lay- 
ers will be separated by a distance equal to the cloud size 
divided by the number of layers minus one, and the de- 
crease of the flux with distance is taken into account from 
layer to layer. We also calculate the flux decrease due to H I 
absorption at Lya, Ly/3, Ly7, LyJ and Lye, for (excitation) 
transitions that are near these wavelengths. The resulting 
population of the ion levels due to photo-exc itati on is calcu- 
lated as described in Vreeswijk et al. (2007) and Vreeswijk 



majority of ionizations are caused by lower energy photons 
(from the ionization threshold of 16.2 eV up to a few hun- 
dred electron-volts). However, the higher-energy photons 
become more important with increasing H I column density, 
as the hydrogen atoms become more and more effective at 
shielding Fe II and other ions from the lower-energy ionizing 
photons. We stress that although we are mentioning Fell 
and Fein here, the same applies to any other ion that is 
included, such as Si II, C n and Cr II. 

Since our calculations involve ionization of Fe II as well 
as the population of excited levels of Fein, we need to take 
into account that after ionization of Fell to Fein, the lat- 
ter will not necessarily be in its ground state immediately. 
Since these are non-standard calculations and rather impor- 
tant for our modelling of the GRB 080310 column densities, 
they are described in Sect. |3.2.2| 

After a packet of photons (AtF v ) has travelled through 
all the layers in the cloud along the sightline, and ionized 
and excited ions along its path, the remaining column den- 
sities (ground state and - if relevant - excited levels) in 
each layer are updated to be used for the next packet com- 
ing through. At each time step, the total ground-state and 
excited-level column densities are determined by simply 
adding up all layer column densities, which can be com- 
pared to measurements. 



et al. (2011) (see also Sect. |3.ip for "each relevant sublayer 
For ionization, the optical depth is: r = Na, and we 
calculate the sub-layer ionization rate as follows: 
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where ov, S ubsheii is the cross section (in cm 2 ) of a partic- 
ular subshell of an ion. The cross sections of the inner shells 



are ta ken from |Verner et al. ( 1993[ ) and Verner & Yakovlev 
( 1995 ) R w hile the cross sectio ns for the outer shells are 



taken from Verner et al. (11996 h . The integration interval, 
set by E t h and F max , is also defined in these references. The 
total ionization rate of each sublayer is obtained by sum- 
ming the rates of the different ion subshclls. The decrease 
in the sublayer column density due to ionization is deter- 
mined USing ^sublayer = -^sublayer ^sublayer ■ For ionS with 

excited-level populations (e.g. Fen and Fein), we assume 
that the ionization cross section is independent of the exci- 
tation level of the ion, i.e. the fraction of ionizations from 
each level is considered to be the same. 

When Fen is ionized by removal of an inner-shell elec- 
tron, i.e. one that is not in the outer 3d or 4s shells, it is 
generally not converted to Fe III. The vacancy created in the 
inner shell is filled by a cascade of radiative (fluorescent) 
and non-radiative (Auger) transitions, which can cause ad- 
ditional electrons to be expelled. The total number of elec- 
trons that are removed following ionization of an electron 
in a particular shell has been calculated for v arious ions and 
ionization states by Kaastra & Mewe (1993), whose results 
we apply in our program. In practice, practically all Fell 
ions that are converted to Fe III are due to outer shell (3c? or 
4s) ionizations. We also note that although higher-energy 
photons (up to several tens of keV) can ionize Fen, the 



9 included in the VizieR data base, at 
Istrasbg.frl with catalog ID: J/A+ AS/109/125 
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